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> ■ Abstract 

X 
b . 

We review studies of an evolution operator C for a discrete Langevin equa- 
tion with a strongly hyperbolic classical dynamics and a Gaussian noise. The 
leading eigenvalue of C yields a physically measurable property of the dynam- 
ical system, the escape rate from the repeller. The spectrum of the evolution 
operator C in the weak noise limit can be computed in several ways. A 
method using a local matrix representation of the operator allows to push 
the corrections to the escape rate up to order eight in the noise expansion 
parameter. These corrections then appear to form a divergent series. Actu- 
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ally, via a cumulant expansion, they relate to analogous divergent series for 
other quantities, the traces of the evolution operators C n . Using an integral 
representation of the evolution operator C, we then investigate the high order 
corrections to the latter traces. Their asymptotic behavior is found to be con- 
trolled by sub-dominant saddle points previously neglected in the perturbative 
expansion, and to be ultimately described by a kind of trace formula. 



I. INTRODUCTION 

In the statistical theory of dynamical systems the development of the densities of particles 
is governed by a corresponding evolution operator. For a repeiler, the leading eigenvalue 
of this operator £ yields a physically measurable property of the dynamical system, the 
escape rate from the repeiler. In the case of deterministic flows, the periodic orbit theory 
]T|fJ] yields explicit and numerically efficient formulas for the spectrum of L as zeros of its 
spectral determinant ||. 

Upon all dynamical evolutions in nature, stochastic processes of various strength have 
an influence. In a series of papers 0-0 the effects of noise on measurable properties such as 
dynamical averages in classical chaotic dynamical systems were systematically accounted. 
The theory developed is closely related to the semi-classical h expansions |8|-[T0| based on 



Gutzwiller's formula for the trace in terms of classical periodic orbits ||11|| , in that both 
are perturbative theories in the noise strength or h, derived from saddle-point expansions 
of a path integral containing a dense set of unstable stationary points. The analogy with 
quantum mechanics and field theory was made explicit in [|| where Feynman diagrams were 
used to find the lowest nontrivial noise corrections to the escape rate. 

An elegant method, inspired by the classical perturbation theory of celestial mechanics, 
was that of smooth conjugations ||. In this approach the neighborhood of each saddle 
point was flattened by an appropriate coordinate transformation, so the focus shifted from 
the original dynamics to the properties of the transformations involved. The expressions 
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obtained for perturbative corrections in this approach were much simpler than those found 
from the equivalent Feynman diagrams. Using these techniques, we were able to extend the 
stochastic perturbation theory to the fourth order in the noise strength. 

In || we developed a third approach, based on construction of an explicit matrix rep- 
resentation of the stochastic evolution operator. The numerical implementation required a 
truncation to finite dimensional matrices, and was less elegant than the smooth conjuga- 
tion method, but made it possible to reach up to order eight in expansion orders. As with 
the previous formulations, it retained the periodic orbit structure, thus inheriting valuable 
information about the dynamics. 

The corrections to the escape rate appeared to be a divergent series in the noise expansion 
parameter. They actually followed, via the cumulant expansion, from the traces of the 
evolution operator £ ||, for which the noise correction series looked similarly divergent. In 
the high order corrections were then worked out analytically for the special case of the 
first trace, Tr(£), as contour integrals asymptotically evaluated by the method of steepest 
descent. Here we further confirm the divergent nature of the high order noise corrections 
for Tr£ ra ; their high order terms are also controlled by sub-dominant saddles, which can be 
interpreted as generalized periodic orbits of some associated discrete Newtonian equations 
of motion. 

In the following sections, first we define the stochastic dynamics, the evolution operator 
and its spectrum. Since this work was inspired by the results obtained for the eigenvalue 
corrections via the matrix representation method outlined in || , for the sake of completeness 
we next show how to obtain a matrix representation of the evolution operator as an expansion 
in terms of the noise strength a, and how to calculate the spectrum of the operator. Then 
we turn towards investigating the behavior of the high order noise corrections to the traces 
of C n , which are responsible for the divergent behavior seen at late terms of the eigenvalue 
corrections. Our key result is fl59|) where the high order noise corrections are converted into 
a trace formula. We give as a numerical example the quartic map considered in 
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II. THE STOCHASTIC EVOLUTION OPERATOR 



In this section we introduce the noisy repeller and its evolution operator. An individual 
trajectory in presence of additive noise is generated by iterating 

Xn+l = f{x n ) +&£n (1) 

where f(x) is a map, £ n a random variable with the normalized distribution P(£), and a 
parameterizes the noise strength. In what follows we shall assume that the mapping f(x) is 
one-dimensional and expanding, and that the £ n are uncorrelated. A density of trajectories 
4>{x) evolves with time on the average as 



4>n+i{y) = {Co <f> n ) (y) = J dx C(y, x)<f> n (x) (2) 

where the evolution operator £ has the general form 

£(x,y) = 5 a (y - f(x)), (3) 

8 a (x) =Js{x- at)P{Qdt = ip (^) . (4) 

For the calculations in this paper, Gaussian weak noise is assumed. The evolution operator 

£( x iy) = 1= e 2a ( 5 ) 

V Z7T<7 



can be expanded, in the perturbative limit a — » 0, as 

oo 

£(x,y) ~ J2 a 2N a 2N 5^ 2N \y - f(x)), (6) 

N=0 

where S^ 2 ^ denotes the 2A-th derivative of the delta distribution and a 2 N = (2 7V A!) _1 (or 
a 2 N = f° r general noise with finite ra-th moment m n ). The map used for concrete 

calculations is the same as in our previous papers, a quartic map on the (0, 1) interval given 
by 



f(x) = 20 
which is also shown on Fig. ([!]). 



I /l 

x 
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(7) 



Throughout the theory developed in previous works j§-0], the periodic orbits of the 
system played a major role. A periodic orbit of length n was defined simply by 

Xj+i = f(x j ), j = l,...,n (8) 
x n+ i = x x . (9) 

(This subscript j will quite generally be defined mod n.) 

For a repeller the leading eigenvalue of the evolution operator yields a physically mea- 
surable property of the dynamical system, the escape rate from the repeller. In the case of 
deterministic flows, the periodic orbit theory yields explicit formulas for the spectrum of £ 
as zeros of its spectral determinant 0]. Our goal here is to explore the dependence of the 
eigenvalues v of £ on the noise strength parameter a. 

The eigenvalues are determined by the eigenvalue condition 

F(cr, v{a)) = det(l - £/v{&)) = (10) 

where F(a,l/z) = det(l — z£) is the spectral determinant of the evolution operator £. 
Computation of such determinants starts with evaluation of the traces of powers of the 
evolution operator 

zC 00 

Tr-—- = Y / C n z n , C n = Tr£ n , (11) 

1 ZL - n=l 

which are then used to compute the cumulants Q n = Q n (C) in the cumulant expansion 

oo 

det(l-z£) = l-Y,QnZ n , (12) 

71=1 

by means of the recursion formula 

Qn — — {C n — C n -iQi — ■ ■ -CiQn-l) (13) 
n 

which follows from the identity 

/ oo n \ 

det(l - z£) = exp - ^ — Tr£" . (14) 



In the next section we show how to compute the cumulants Q n using a local matrix repre- 
sentation of the evolution operator. 



III. THE SPECTRUM OF THE EVOLUTION OPERATOR 



As the mapping f(x) is expanding by assumption, the evolution operator (fj) smoothes 
the initial distribution <p(x). Hence it is natural to assume that the distribution <p n (x) is 
analytic, and represent it as a Taylor series, intuition being that the action of C will smooth 
out fine detail in initial distributions and the expansion of (j) n {x) will be dominated by the 
leading terms in the series. 

An analytic function g(x) has a Taylor series expansion 



oo g,m Qm 

q(x) = — - — — q(y) 



y=0 



Following H.H. Rugh [12| we now define the matrix (m, m' = 0, 1, 2, ...) 

L is a matrix representation of C which maps the x m component of the density of trajectories 
4> n (x) in (0) to the y m ' component of the density 4> n+ i(y), with y = f(x). The desired traces 
can now be evaluated as traces of the matrix representation L, Tr£ n = TrL™ . As L is 
infinite dimensional, in actual computations we have to truncate it to a given finite order. 
For expanding flows the structure of L is such that its finite truncations give very accurate 
spectra. 

Our next task is to evaluate the matrix elements of L. Traces of powers of the evolution 
operator C n are now also a power series in a, with contributions composed of 8^ n \f{x a ) — 
x a+ \) segments. The contribution is non-vanishing only if the sequence Xi,X2, ...,x n ,x n+ i = 
x\ is a periodic orbit of the deterministic map f(x). Thus the series expansion of Tr£ n has 
support on all periodic points x a = x a+n of period n, f n (x a ) = x a ; the skeleton of periodic 
points of the deterministic problem also serves to describe the weakly stochastic flows. The 
contribution of the x a neighborhood is best presented by introducing a coordinate system 
(p a centered on the cycle points, and the operator ([|) centered on the a-th cycle point 

x a -^Xa + <Pa, a = l,...,n p (16) 



/a(0) =f(x a + <f>) (17) 
C a ((f) a+ i, (f) a ) = C(x a+ i + (f>a+l,X a + a ) . (18) 

The weak noise expansion for the a-th segment operator is given by 

oo 

C(0',0) = ^(-a) 2Ar a 2iV( 5( 27V )(0 / + x a+1 -/ a (0)). (19) 

JV=0 

Repeating the steps that led to (p~5|) we construct the local matrix representation of C a 
centered on the x a x a+ i segment of the deterministic trajectory 



m 



Qm 1 

'La) , = — — t / d<b £J<b' <b)^-r 



ml 



(20) 

0'=O 



= E (~C r ) 2N Ct2N(Ba)m'+N,m ■ (21) 

N=max(m—m' ,0) 

Due to its simple dependence on the Dirac delta function, Ba can expressed in terms of 
derivatives of the inverse of f a (<f>)'- 

gr, 



[Ba) nm = J d(f)5((f)' + x a+1 - fa{4>))— 
d n {f-\x a+1 + cj>')-x a y 



(22) 

4>'=o 



d<P ln m\\f' a {f^{x a+l +4>'))\ 
sign(^) d n+ \F a {<PT l+1 



(23) 

(f>'=a 



(24) 

<t>'=o 



(m+l)\ d<p' n+1 

where we introduced the shorthand notation J r a (4 > ') — fa l ( x a+i + 4 ) ') ~ x a- The matrix 
elements can be easily worked out explicitly using fl24j) . In we show that Ba is a lower 
triangular matrix, in which the diagonal terms drop off exponentially, and the terms below 
the diagonal fall off even faster. This way truncating Ba is justified, as truncating the matrix 
to a finite one introduces only exponentially small errors. 

In the local matrix approximation the traces of evolution operator are approximated by 

oo oo 

Tr£"L ddlcs = E np E WW = E C n , N a 2N , (25) 

P r=l N=0 

where TrLp = TrL„ p Li2 • • ■ Li is the contribution of the p cycle, and the power series in a 2N 
follows from the expansion (PH ) of La in terms of Ba. The traces of L n evaluated by fl2"ip 
yield a series in a 2N , and the a 2N coefficients Q n ,N in the cumulant expansion 



F = det(l - zC) = 1 - £ Qn,NZ n a 2N (26) 

n=l Af=0 

are then obtained recursively from the traces, as in ([13]): 

i / n-l iV \ 

Qn,V = — ( Cn.TV — X/ X! Qk,N-lC n -k,l J • (27) 
71 V fc=l Z=0 / 

From the cumulants, by manipulating formal Taylor series, we can calculate perturbative 
corrections to the eigenstates. 

Figure (0) shows the cumulants obtained in the numerical tests and table (|) shows 
the eigenvalue corrections computed from the cumulants. Both the cumulants and the 
eigenvalue corrections exhibit a super-exponential convergence with the truncation cycle 
length n, whereas the corrections appear to form a divergent series in the noise parameter 
a. 



IV. TRACE FORMULA FOR NOISE CORRECTIONS 

With as many as eight orders of perturbation theory, we should now turn towards in- 
vestigating the asymptotic nature of these perturbative expansions. The corrections to the 
escape rate are calculated via the cumulant expansion from the traces of the evolution op- 
erators C n . These traces are responsible for the asymptotic behavior since they themselves 
apparently have divergent high order noise corrections. Thus in this section we investigate 
the behavior of the late terms in the noise expansion series of the traces Tr(£ n ). With the 
help of an integral representation of the operator we shall be able to transform Tr(£ n ) to 
contour integrals, and then we will evaluate these integrals in the saddle point approximation 
to arrive to a formula analogous to the usual trace formulas arising in quantum chaos. 

The trace of C n can be expressed as 

1 r s(s) 

Tr£ n = — == / dxidx2---dx n e , (28) 

(V27T0-) n J - 

where 



S 



1 n 

S(x) = -J2(x j+1 -f(x j ))\ (29) 

Z 3 = 1 

x = (xi, • ■ ■ ,x n ), with x n+ i = xi. (30) 

In order to give a deeper insight on the forthcoming calculations, we draw a correspon- 
dence between our system and a discrete Hamiltonian mechanics, with the S defined above 
playing the role of the classical action. According to (p9]), the least action principle requires 

X j ~ f( X j-l) - f'(Xj)(x j+1 - f(Xj)) = 0. (31) 

We define 

P= (Pi,' ■ ■ ,Pn), Pj ■= Xj - f(Xj-i), (32) 

the quantity corresponding to the momentum in the classical mechanics. From (|3~ID we 
obtain 

xj+i = f(xj) +Pj+i, (33) 

^ ^ jkr (34) 

which are the equations corresponding to the classical Newtonian equations of motion. The 
generalized periodic orbits of length n are those orbits, which obey these equations and 
x n+ i = X\, p n +i = p n . Those generalized periodic orbits which have non-zero momentum will 
control the asymptotic behavior of the corrections to Tr£ n as we shall demonstrate later. The 
original periodic orbits defined by are those with zero momentum. The generalized 

periodic orbits with non-zero momentum and the original periodic orbits proliferate with 
growing n as suggested by Fig [| 

We introduce an integral representation of the noisy kernel, which will be of great use in 
the later calculations: 

1 (y-ft*)) 2 

£{x,y) = e = 

V Zncr 



— I dke-^ +ik{ -y-^\ (35) 
2n J 



Using this new integral representation, 



Tr£ n 
1 

(2n) n 



dk n dx n e-^U k > +i ^U k ^+i-K x ^. 



or equivalently 



Tr£ n 
1 

(2ir) n 



dk n / dp n J n (p) e' 



where J n (p) denotes the Jacobian D(x)/D(p). Since 

' r dk n e^U k ^ = f[S{ Pj ) 

3=1 



(27T) 



we can reduce (]37f) to 



Tr£ n = J dp n J n (p) e 



— iin 



n <%: 



Pi=0 



with A n denoting the Laplacian 

A r , = 



dp\ dpi "' dp 2 n 



We focus on the Taylor expansion of fl55|) in the noise parameter: 



Tr£ n = (Trrj^fj^ 

^=0 



(Tr£ n ; 



1 (A r 



N 



' J'n 



Pj =0 



2 N N] "» 

The iV-th power of the Laplacian in the equation above can be written as 

31 



(A r 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 
(42) 



(43) 



where Sji is the Kronecker-delta (making the sum actually finite). With the help of the 
multidimensional residue formula from complex calculus |13| 
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Qm + ...+n k j fg\ 

dzT...dz n k h 
ni\...n k \ 

(2iri) k 



we obtain 



(Tr£ r 



(&-*i) ni+1 ...(&-**) nfc+1 ' 

(2ji)! (2 3n ))\ 



(44) 



1 



E 



x (2ni)»2» ilr ^ n=0 ./,! "' in! 

</n (P) dpi...dp r , 



2ji+l 2j„+l • 



(45) 



These contours encircle the pj = points positively. The integrals can be transformed back 
to contour integrals in the original Xj variables, and the contours will be placed around the 
original periodic orbits of the system defined by (|8|-|9|), since it is these orbits which fulfill 
the pj = conditions. 



(Tr/7 



1 



E 



(2ji)< (2j n ))\ 



(27ri)»2"^ n=0 ./.! jj 



xS 



Jk 



dxi...dx n 



{xi - f{x n )fn+\..( Xn - f(x n ^))^+y (46) 
From now on we shall restrict our calculations to the asymptotic large-iV limit. We will 
replace the summations in ( fj5|) by integrals and then use the saddle-point method to get a 
compact formula for (Tr£ ra )jv. Then we may consistently approximate the factorials via the 
Stirling formula [TM as 



(2j k y. [^y^v^Tk 



3k\ 



22jk+l/2jjk e -jk 



2l/2 e 2(ln2)j fe +i fc lni fc -j fe _ 



(47) 



Using (|47| ) and an integral representation of the delta function we get 

9f-JV oo 

(TrC n ) N ~ t r Y dti ..A dx x ...dx n 

n n n 

x exp it(N - ]T Jk ) + (2 In 2 - 1) £ j fc + £ j fc In Jfc 



L fc=l fe=l 

n 

^ln(x fc -/(x fc _ 1 ))(2 Jfc + l) 



fe=i 



(48) 



Now we replace j k with the new variables y k = % and in the asymptotic (N large) limit 
approximate the summations over j k with integrals over y k as 



(Tr£ r 



2 %-N N n fOQ 



>N 



, , dy x ... dy n dti ... f dx 1 ...dx r , 

{27it) n 27T Jo Jo J Jci Jc n 



x exp 



it{N -NY J Vk)+ N{2 In 2 - 1) £ y k 

k=l k=l 



+N Vk HNy k ) + E ln (^ - f(?k-i))(2Ny k + 1) 

k=l k=l 



We evaluate the y integrals with the saddle-point method to get 

2 



(49) 



(Tr£ n ; 



N 



(2vr)2i n 27r 



cxp 



it L/V + 



dt 


fork 






n\ 


i,S] 


2j 


~ e 2. 



dx i • • • dx yi 



(50) 



Next we implement the saddle-point method to the integral in t as well, asymptotically 
resulting in 

ivV (2N)\ 



(Tr£ n ) 



iV 



2 2iV+ 5(27r) Ii * i ^+ 1 N\ 
x [dx n e -( N+ ^ s K 



(51) 



The last step is to evaluate the contour integrals in the x k variables. We deform the contours 
until the saddle-points are reached and the contours run along the paths of steepest descent. 
The leading contribution comes from the saddle-points, which fulfill the following equation 

1 



S 



■>■) - ./•(•'-.; ,) - (3+1 - /(**))/'(**) 



0. 



(52) 



By comparing ([52]) and (|3T| ) one can see that the saddle-points are all generalized periodic 
orbits of the system. Since the contours ran originally around the orbits with zero mo- 
mentum, these do not come into account as saddle-points. The second derivative matrix 



is 



12 



(53) 

where D 2 S denotes the second derivative matrix of S 

^ = (54) 

This would be the matrix to deal with if we were taking the saddle-point approximation of 



directly. We reorganize the prefactor in (51) with the use of the Stirling formula [14 



then the result of the saddle-point integration comes out as 



N ^ 



(Tr£ n ; 

E gV r(AT + i) s;» 
tt 2 " (iV + f) T s /dtiD% 

For n — 1 this formula restores the result of [0] as it should. 

Finally we draw attention to the close connection between the generalized periodic orbits 
of the system and D 2 S. The stability matrix of a generalized periodic orbit p is expressed 
as 

M p = M{n) ■ M{n - 1) • M(n - 2) ■ • ■ M(l) (56) 

M{k) = (fW ~ W$frf" rk)\ (57) 

The determinant of D 2 S can be expressed with the help of the stability matrix as 

det D 2 S P = det(M p - 1). (58) 
This way we reformulate fl5"5|) as 

N 1 ^ 1 T(N + i) „ e ~NlnS p 

- 2n A E 7= 7TT ' ( 59 ) 



iV+f ) 2 g.p.o. ^/det(l-M p ) 



where the summation runs over generalized periodic orbits, with non-zero momentum. This 
expression, fully analogous to a trace formula, is our main result. 

Finally, we provide a numerical test of our high order estimates. In [[7] we developed 
a contour integral method to calculate high order noise corrections to the trace of £, and 
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obtained a very good agreement between the exact results and a formula which is just the 
approximation of fl5"5| ) in the n — 1 case. Now we can similarly verify the high order noise 
corrections to the trace of C 2 . Fig. |] plots, as a function of N, the ratios of the approx- 
imations for (Tt£)n and (TtC 2 )n obtained from ([55]) to the exact results ([46]) computed 
using direct numerical evaluations of the contour integrals therein. (For higher powers of 
C, the proliferation of the orbits forces the contours to run close to singular points, making 
efficiently precise numerical integrations extremely time-consuming). 

In summary, we have studied the evolution operator for a discrete Langevin equation 
with a strongly hyperbolic classical dynamics and a Gaussian noise distribution. Motivated 
by the divergent growth of the high order eigenvalue corrections of the evolution operator, 
we investigated the behavior of the high order noise corrections to the traces of C n . Using an 
integral representation of the evolution operator C, we found that the asymptotic behavior 
of the corrections to the trace of C n is governed by sub-dominant quantities previously 
neglected in the perturbative expansion, and a full-fledged trace formula can be derived for 
the late terms in the noise expansion series of the trace of L n . 

G.V. and G.P. gratefully acknowledge the financial support of the Hungarian Ministry 
of Education, EC Human Potential Programme , OTKA T25866. G.P., G.V. and A.V. were 
also partially supported by the French Ministere des Affaires Etrangeres. 
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FIGURES 

FIG. 1. The map (7) on the [0,1] interval. 

FIG. 2. The perturbative corrections to the cumulants Q„, : n plotted as a function of cycle 
length n (for perturbation orders N = 0,2,4,6,8); all exhibit super-exponential convergence in n. 

FIG. 3. The shortest original and generalized periodic orbits of the map (7). Squares mark 
original periodic points, dots mark generalized periodic points. Large symbols indicate orbits of 
length one, small symbols indicate orbits of length two. 

FIG. 4. The ratios of (Tr£)jv and (Tr£ 2 )7v as calculated via the asymptotic formula (55) to 
their values computed by numerical integration of formula (46) . 
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FIG.2. 
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FIG.4. 
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TABLES 



n u 



V2 



1 0.308 

2 0.37140 

3 0.3711096 



0.42 

1.422 

1.43555 



4 0.371110995255 1.435811262 



2.2 
32.97 
36.326 
36.3583777 



17.4 
1573.3 
2072.9 
2076.479 



168.0 
112699.9 
189029.0 
189298.8 



5 0.371110995234863 1.43581124819737 36.35837123374 2076.4770492 189298.12802 

6 0.371110995234863 1.43581124819749 36.358371233836 2076.47704933320 189298.128042526 
TABLE I. Significant digits of the leading deterministic eigenvalue vq, and of its a 2 ,---,a s 

perturbative coefficients zv 2 , • • • , f 8 , calculated from the cumulant expansion of the spectral deter- 
minant as a function of the cycle truncation length n. Note the super-exponential convergence of 
all coefficients. 
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